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Abstract 

X> ■ 

, The gluing technique is used to construct hypersurfaces in Euclidean space having approx- 

pLn ' imately constant prescribed mean curvature. These surfaces are perturbations of unions of 

( ^-> ) \ finitely many spheres of the same radius assembled end-to-end along a line segment. The con- 

■ dition on the existence of these hypersurfaces is the vanishing of the sum of certain integral 
moments of the spheres with respect the prescribed mean curvature function. 

Q ■ 1 Introduction 

•i— > . 

In the recent paper [1J, Butscher and Mazzeo have constructed examples of constant mean curvature 
(CMC) hypersurfaces in a Riemannian manifold M with axial symmetry by gluing together small 
spheres positioned end-to-end along a geodesic 7. The examples they have constructed have very 

' large mean curvature 2/r and lie within a distance 0(r) of either a segment or a ray of 7; hence 

we say that these surfaces condense to the appropriate subset of 7. Such surfaces cannot exist 
in Euclidean space, and their existence relies on the fact that the gradient of the ambient scalar 

I curvature of M acts as a 'friction term' which permits the usual analytic gluing construction (akin 

to the classical gluing constructions pioneered by Kapouleas [2], [1]) to be carried out. The purpose 
of this paper is to show the same techniques used in [1] can be adapted in a straightforward manner 

■ to show that a similar construction is possible in a much simpler yet fairly general context: that of 
^ ■ hypersurfaces having near-constant prescribed mean curvature in Euclidean space. The essence of 

the gluing construction carried out herein therefore lies in identifying and appropriately exploiting 
the analogous 'friction term' appearing in this setting. 

Let F : M n+1 x TM n+1 ->Rbea given, fixed smooth function. For simplicity and to maintain 
the parallel with [I], we will assume that F has cylindrical symmetry in the following sense. Endow 
R n+1 with coordinates (x°, x , . . . , x n ) and let G C 0(n+l) be the set of orthogonal transformations 
that fix the x°-axis. Each rotation R G G acts on TR n+1 via the differential R* : TR n+1 -> TW a+1 . 
We will now demand that F(R{p),R*V p ) = F(p, V p ) for all (p, V p ) G R n+1 x TR n+1 . The prescribed 
mean curvature problem that will be solved in this paper is to find, for every sufficiently small 
r G K+, a G-invariant hypersurface S r which satisfies 

£T[E r ](p) = 2 + r 2 F(p,Nx r (p)) VpGS r (1) 
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where i?[S r ] is the mean curvature of S r and Ns r is the unit normal vector field of S r . Furthermore, 
this hypersurface will be constructed by gluing together a finite number K of spheres of radius one 
(and thus of mean curvature exactly equal to two) whose centres lie on the x°-axis using small 
catenoidal necks having the x°-axis as their axes of symmetry. 

In order to properly state the Main Theorem, we must make the following definition, which is 
meant to capture the most important effect of the prescribed mean curvature function F on the 
surface that we're attempting to construct in this paper. 

Definition 1. Let S be a compact surface in W l+ . The F -moment of S is the quantity 

MS) ■= [ F(x,N s (x))JdVol s 
Js 

where N$ is the unit normal vector field of S and dVol^ is the induced volume form of S, while 
J : S -> M is defined by J(x) := N s (x)} for x £ S. 

Now let p®(s) := (s + 2(k — 1),0, ... ,0) and consider the spheres Sk(s) := dBi(p^.(s)). These 
spheres are positioned along the x°-axis in such a way that each Sf~(s) makes tangential contact 
with Sfc±i(s). The following theorem will be proved in this paper. 

Main Theorem. Suppose that there is sq € M. such that 

• the F -moments of the spheres Sk{so) satisfy ^2k=i M-F^fc^o)) = 

• the function s i— ► Ylk=i ^F^ki 8 )) has non-vanishing derivative at s = sq, 

then there for all sufficiently small r > 0, there is a smooth, embedded hypersurface E r which is a 
small perturbation o/ljfLi S^sq) that satisfies the prescribed mean curvature equation ([1]). 

It is easy to find a situation in which the conditions of the Main Theorem hold. For example: 
if F(-, •) is such that iJ 1 p(dBi(x°,x 1 , . . . ,x n )) is negative whenever x° is sufficiently negative and 
positive whenever x° is sufficiently positive, then the mean value theorem asserts that a zero of the 
function s i— > Ylk=x I^FiSkis)) can be found. And if also F(x, •) is monotone as a function of x°, 
then this function will have non-zero derivative. 

An application of the Main Theorem, and indeed an inspiration for it, is the earlier work by 
Kapouleas on slowly rotating assemblies of water droplets [3]. In this case, the prescribed mean 
curvature function F : l n+1 x TIR n+1 -> R takes the form F{p,N Tlr {p)) := C(uj)(p ) 2 where 
p := (p ,^ 1 , . . .p n ) and C(u>) depends on the angular velocity us. The prescribed mean curvature 
equation now approximates the effect of centrifugal force on the surface S r when u is small. One 
of the assemblies of water droplets that Kapouleas constructs is exactly as described in the Main 
Theorem. (He constructs many other, more complex, and less symmetrical assemblies as well.) 

Another application of the Main Theorem is for understanding the possible shapes an electrically 
charged soap film can adopt in the presence of a weak, axially symmetric electric field. In this case, 
the equation satisfied by the surface adopted by the soap film is exactly ([!]), where the prescribed 
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mean curvature function F : I n+1 xTK" +1 — > R takes the form F(p, Ns r (p)) := -C(V0(p), iV Er (p)} 
while : R™ +1 — > R is the electric potential and C is a constant. We can see why this is so by 
writing the total energy of the soap film as the sum of a surface area term and a term proportional 
to the surface integral of eft, and then computing the Euler-Lagrange equation for the variation of 
this energy subject to the constraint that the volume enclosed by the surface remains constant. 
If we now assume that <p is such that the existence conditions of the Main Theorem hold, then 
the Main Theorem asserts that K spherical, electrically charged soap films connected by small 
catenoidal necks can be held in equilibrium at special points in space by the electric field. 

2 The Approximate Solution 

To construct an approximate solution for the Main Theorem, we use essentially exactly the same 
procedure as in [TJ §3.1]. This will be outlined here very briefly for the convenience of the reader. 
The presentation is given for the dimension n = 2 for simplicity; everything that follows can be 
easily adapted to the (n + l)-dimensional setting. 

Endow R n+1 with coordinates x , . . . , x n ) and let 7 be the arc-length parametrization of 
the x°-axis with 7(0) = (0, 0, ...,0). We will construct an approximate solution for the Main 
Theorem out of K spheres of radius one as follows. Choose a localization parameter s € R and 
small separation parameters a%, . . . , ctk-i £ R+- Define s\ := s and Sk := s + 2(k — 1) + Y2i=i a l 
for k = 2, ... K and set pk ■= 7(sfe) and p^ := 7(5^ ± 1). Define the spheres Sk ■= dB\{pk). These 
spheres will now be joined together according to the following three steps. 

Step 1. The first step is to replace each Sk with the surface Sk obtained by taking the normal 
graph of a specially chosen function Gk over Sk \ [B Pk (p~^) U B Pk (p7)] where pk € (0, 1) is a small 
radius as yet to be determined. The functions we use for this purpose satisfy the equations 

• C{G k ) = e+5(p+) + £76{p7 ) + A k Jk if k = 2, ... K - 1 

• C(d) = ef5(j>+) + A X J X if k = 1 

• C{G K ) = £k 6 (Pk) + A k j k if fc = K 

where C := A g 2 + 2 is the linearized mean curvature operator of the unit sphere, the small scale 
parameters are yet to be determined and 5(q) is the Dirac ^-function centered at q, while 
Jk := (-7^5, Ng k ) is the sole G- invariant function in the kernel of C normalized to have unit L 2 - 
norm, and Ak is chosen to ensure L 2 -orthogonality to Jk- Of course Jk = x °\g k > t ne restriction of 
the x° coordinate function to Sk- 

Step 2. Let W be the catenoid, i.e. the unique complete minimal surface of revolution whose axis 
of symmetry is 7 and whose waist lies in the (x 1 , x 2 )-plane. The next step is to find the truncated 
and re-scaled catenoidal neck of the form Wk ■= B p > (p b k ) n [e^ W +p\ + (Sk, 0, 0)] that fits optimally 
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in the space between Sk and Sk+i for k = 2, . . . , K — 1. Here > is a small scale parameter and 
is a point between p k and p^ +1 that are determined by the optimal fitting procedure while 6k 
is a small displacement parameter that takes Wk away from its optimal location and p' k is a small 
radius as yet to be determined. The optimal fit is obtained by matching the asymptotic expansions 
of the functions giving Sk D B p ^ (p k ) and Sk+i H -B p j_ (p^,) and Wfc as graphs over the translate of the 
(x 1 , x 2 )-plane passing through p, exactly as in [Tj §3.1]. One particularly important outcome of the 
matching is that Sk from the previous step, as well as ef and p\ are all uniquely determined by ak- 
in fact, an invertible relationship of the form cr^ := Afc(e^) holds, with Ak(£k) = 0(ek\ log(efc)|). 
Finally, we find that we must choose Pk-,p'k = ^( e ^ 4 ) to ensure the optimal fit between the necks 
and the perturbed spheres. 

Step 3. The final step is to use cut-off functions to smoothly glue the neck Wk{&k) into the space 
between Sk and §k+i- The interpolating region is the annulus B p i (p k ) \ B p ^ 2 (Pk)- 

In this way we obtain a family of surfaces depending on the a, 6 and s parameters. 

Definition 2. Let K be given. The approximate solution with parameters a := {<n, . . . ,o~k—i} 
and 6 := {6%, . . . , 6k-i} and s is the surface given by 





" K 




'K-l 


S r (cj, 6, s) := 




U 


U *k 




.fc=l 




.k=l 



3 Solving the Projected Problem 

We now proceed to solve the equation ([TJ up to a finite-dimensional error term by perturbing the 
approximate solution constructed in the previous section. The required analysis is in most respects 
identical to or less involved than the analysis found in [J, §4 - §6] and will thus again only be 
abbreviated here for the sake of the reader. The outcome will be a surface Tj\{a : 6, s) satisfying 
i?[x|(<7, 6, s)] — 2 — r2 ( a 8 s) ^ wnere VV is a finite-dimensional space of functions that will 
be defined precisely below. It arises because the linearized mean curvature operator, which governs 
the solvability of ([2]) , possesses a finite-dimensional approximate kernel consisting of eigenfunctions 
corresponding to small eigenvalues. These small eigenvalues make it impossible to implement a 
convergent algorithm for prescribing the components of the mean curvature of the approximate 
solution lying in W. 

Function spaces. We first define the weighted Holder spaces in which the analysis will be carried 
out. These are essentially the same weighted spaces as in [Tj §4], namely the spaces C^' Q (E r (c7, 6, s)) 
consisting of all functions on S r (<7, 6, s) where the rate of growth in the neck regions of 

X r (<7, 6, s) is controlled by the parameter v. Choose some fixed, small < ft < 1 and define a 
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weight function £ r : X r (c, 5, s) — > R as 



\\x\\ p = (x , x) £ B R j 2 {Pk) f° r some k 

Interpolation p E Bn{p^ k ) \ B^j 2 {p\) f° r some k 
1 elsewhere 



where the interpolation is such that £ r is smooth and monotone in the region of interpolation, has 
appropriately bounded derivatives, and is G-invariant. Now for any open set U C E r (<r, 5, s), define 

k 
i=0 

where | • \ou is the supremum norm on IA and Wot U is the a-H61der coefficient on U. This is the 
norm that will be used in the C v ,a (jb r {o~, 5, s)) spaces. 

The equation to solve. Let [i : G^' a ($] r ((7, 5, s)) — > Emb(f^ r (a, 5, s), R" +1 ) be the exponential 
map of S r (<7, 5, s) in the direction of the unit normal vector field of S r (cr, 5, s). Hence /x^(S r (cr, 5, s)) 
is the normal deformation of S r (cr, 5, s) generated by / G C^' a (S r (cr, 5, s)). The equation 

H[^ f (t r (a,5,s))] = 2 + r 2 F o x N f{tMs)) ) (2) 



selects / € c3' Q (£ r (<7, 5, s)) so that /i/(S r (cr, 5, s)) satisfies equation ([T]). In addition, the function 
/ will be assumed G-invariant. Define the operator 

$ r , w : C 2 u > a (t r (a,5,s)) -> C°l Q 2 (S r ((7,«5, S )) 
*r,.,«T,« (/) := H [M/(Sr (o-, <5, a))] - 2 - r 2 ^(/) 

where :=Fo( 

Z 2 / x ^(s r (u,5,s)) J- ^he linearization of $V )0 - )( 5 )S at zero is given by 
C := Ddv As (0) = A + ||£|| 2 + r^D^Guo, A^ M , s) ) • /^(..i,.) " ^F^o, %. M , s) ) ■ V/) 

where DiF and D2-F are the derivatives of F in its first and second slots and B := B[E r (o~,5,s)] 
is the second fundamental form of E r (<r, 5, s). On the k th spherical part of S r (cr, 5, s), the operator 
£ is a small perturbation of Ck '■= As fc + 2 which is the linearized mean curvature operator of the 
sphere Sk- Let once again be the G-invariant function in its kernel and define the space 

VV := sp&n{xext,kJk , Xext,k£k(Vk) ■ k = 1, . . . , K - 1} U {Xezt,K<//<} • 

Here Xexi,fc is a smooth cut-off function supported on Sk and is a smooth cut-off function 
supported on the transition region between the k th neck and with the property that the support 
of Vrjk and Vxext,k do not overlap. 

We now prove the following theorem. Let e := max{ei, . . . , sk-i} and 5 := max{(5i, . . . , 5k -1} 
and we will assume that e = 0(r 2 ) and 5 = O(r), which will be justified a posteriori. 
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Theorem 3. If r is sufficiently small, then there exists f := f r (a,5,s) G Ci/ a (I] r (<7, S, s)) wit/i 
1/ G (1, 2) so i/tai 

*W,.(/)GW. (3) 

The estimate \ f\ c 2, a < Cr 2 ZioWs /or i/ie function f , where the constant C is independent of r. 
Finally, the mapping (o~,5,s) i— > f r (o~,5,s) is smooth in the sense of Banach spaces. 

Proof. As in [I], we will use a fixed-point argument to solve the equation & r ,a,6,s{f) € W for 
a function / G C^ ,Q: (S r (o", 5, s)) with f G (1,2). The fixed-point argument follows from three 
steps: an estimate of the size of $ r ,s,a,S,s(Q)l the construction of a bounded parametrix 1Z satisfying 
C o 1Z = id + £ where £ : C^^^rO 7 , s )) ~* W; an d an estimate of the non-linear part of the 
operator & r ,a,S,s- Each of these steps is given in great detail in pp. Thus all that is needed here is 
to point out how the analysis of [1] applies to the present situation. 

Step 1. We begin with the estimate of | $^,5^(0)1 £,0,0 , which is the amount that the approximate 
solution £ r (o~, S, s) deviates from being an actual solution of equation This is done by adapting 
[TJ Prop. 13]. In fact, by combining Steps 1, 2 and 4 of the estimate of H[T, r (a, S, s)] — 2 in the 
C v —2 f° r v ^ (1) 2) found in [TJ Prop. 13], with a straightforward estimate for the C v ^ 2 norm of the 
r 2 J- term, we find that 

|$r«ri.(0)U« < Cmaxjr 2 , £ 3/2-3^/ 4) fe l-3,/4| < Cr 2 

for some constant C independent of r. 

Step 2. We now find a parametrix 1Z : C°^* 2 (£v(<r, 5, s)) — > Cv' a (£ r (er, s)) satisfying £olZ = 
id + £ where £ : C^* 2 (E r (<r, <J, s)) — > VV. As in [TJ Prop. 15], this is done by first constructing 
an approximate parametrix by patching together parametrices for the linearized mean curvature 
operator of each sphere with parametrices for the linearized mean curvature operator of each neck; 
and then iterating to produce an exact parametrix plus an error term in W in the limit. The 
difference here is that the terms coming from the non-Euclidean background metric in [1, Prop. 15] 
must be replaced by the r 2 T term. The same result holds because this term can easily be shown to 
satisfy the right estimates. In fact, 7Z and £ satisfy the estimate \1Z(w)\ c 2, a + \E(w)\ c a, a < C\w\ c o, a 
for all w G C°lf 2 (^r(0', <^> s ))i where C is a constant independent of r. 

Step 3. We define the quadratic and higher remainder term of the operator & ra 6s as 

Q : C^ a (± r (a, 5, s)) - C°L Q 2 (S r (a, 6, s)) 

Q(f) ■= *r,v,sM) ~ *W,.(0) " OJ) 

The estimates for the Ct ,a norm of Q can be found exactly as in [1, Prop. 18] with the terms 
coming from the non-Euclidean background metric replaced by the r 2 J- term. Then there exists 
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M > so that if fx, / 2 S Cv' a (E r (a, 5, s)) for v £ (1, 2) and satisfying l/i^a + |/2| c 2,a < M, then 
|Q(/i) - Q(/2)| c o L « < C|/i -/ 2 | C 2, a max{|/i| C 2,a,|/2| c v*} 

where C is a constant independent of r. Once again, this works because the r 2 T term can easily 
be shown to satisfy the right estimates. 

Step 4. We can now solve the CMC equation up to a finite-dimensional error term by implement- 
ing a fixed-point argument based on the parametrix constructed in Step 2 as well as the estimates 
we have computed so far. Let E := $ r ,o-,<5,s(0) and use the Ansatz f := lZ(w — E) to convert the 
equation $V,<r,<5,s(/) £ into the fixed point problem w — M rtCT ,6,s(w) G W where 

K,a,sA w ) ■= ~Q ° K{w - E) + r 2 T o TZ(w - E) 
The estimates that have been established up to now give us the estimate 

\N r {wi) — N r {W2)\ o 0,a < Cr 2 \w\ — W2\ c 0,a 

for w belonging to a ball of radius 0(r 2 ) about zero in C^ 2 (S r (cr, 5, s)), where C is independent 
of r. Hence M r is a contraction mapping on this ball if r is sufficiently small, and a solution of 
([3]) satisfying the desired estimate can be found. The smooth dependence of this solution on the 
parameters (a, 5, s) is a consequence of the fixed-point process. □ 



4 Force Balancing Arguments and the Proof of the Main Theorem 



When r is sufficiently small, we have now found a function f r (a,5) £ C* ,a (T> r (a, 5, s)) for each 
(a, 5, s) so that 

H[nf r { a ,S) (S r (o", 5, s))] - 2 - r 2 J"(/ r ((T, 5, s)) = £ r (a, 6, s) 

where £ r (a,5,s) is an error term belonging to the finite-dimensional space W depending on the 
free parameters (<r, 5, s). The corresponding surface that satisfies the prescribed mean curvature 
condition up to finite-dimensional error is Y,l(a,5,s) := fj,j r ^ a ^^(T, r (a, 5, s)). 

To complete the proof of the main theorem, we must show that it is possible to find a value 
of (a, 5, s) for which these error terms vanish identically. As in [H §7.2], we consider the balancing 
map B r : R 2 ^" 1 -> R 2X_1 defined by 



B r (a, 5, s) := (^i (£ r (cr, 5, s)) tt2K-i (£ r (cr, 5, s)) J (4) 
where iT2k+i '■ VV — > M and 7T2/;; : VV — > K are the L 2 projection operators given by 



7T2fc(e) := 



/ „ e-Xneck,kh and vr 2fc+ i(e) := / e • x' ex t,k J k 

JT,l(cr,S,s) JH$.(a,8,s) 
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where x' ex t k ^ s a cut-off function supported on the k th spherical region, Xneck,k is a cut-off function 
supported on the k th neck and transition region, and I k the Jacobi field of the neck coming from 
translation along the neck axis. This is an odd, bounded function with respect to the centre of 
the neck. Note that B r is a smooth map between finite-dimensional vector spaces by virtue of the 
fact that the dependence of the solution f r (o~,5,s) on (a,5,s) is smooth and the mean curvature 
operator is a smooth map of the Banach spaces upon which it is defined. The following lemma 
proves that 7r(e) = implies that e = 0, and its proof is a straightforward computation. (In order to 
make this work out, we must choose the cut-off functions properly: we must have overlap between 
the supports of Xneck,k, Xext,k and Xext,k+i and no overlap between the supports of x' ex t,k and Vk-) 

Lemma 4. Choose e G VV as e = J2k=i a kXext,kJk + J2k=i b kXext,k^k(m) f or h G R. Then 

^2k(e) = Cib k + C^e^afc+i - e 3 k /2 a k ) 
7T2fc+l(e) = C 2 a k 

where Ci,C[,C 2 are independent of r and (o~,5,s). 

We must now show that B r (a,5,s) can be controlled by the initial geometry of Y, r (o~,5, s), at 
least to lowest order in r. The calculations are similar to those found in [lj §7.2] except with 
the contributions from the ambient background geometry replaced by a contribution from the 
prescribed mean curvature in the form of the F-moments of the spheres making up S r (<r, 8, s). 

The highest-order part of £ r (o~, 5, s) involves the F-moments of the spherical constituents S k 
of S r (cr, 5, s) as follows. Set fi k (cr,s) := /J,F(S k ) — this depends on s and ai,...,a k because the 
location of the centre of S k is determined by these parameters. Let us continue to assume that 
= 0(r 2 ) and 5 k = 0(r) for each k. This will be justified shortly. 

Lemma 5. The quantity £ r (o~,5,s) satisfies the formulce 

ir 2k {£ r (a, 5, s)) = + 0{r 2+2v ) (5a) 

C 2 e 1 -r 2 f i 1 (a,s) + 0(r 4 ) k = 

C 2 {e k+l -e k ) -r 2 fi k+1 (a,s) + 0{r 4 ) 0<k<K-l (5b) 

-C 2 e K -r 2 fi K (a,s) + 0(r 4 ) k = K-l 
where Ci,C 2 are constants independent of r,o~,6,s. 



n 2 k+i(£r(o-, 8,s)) 



Proof. Set Sr := Sr(cr, 5, s) and S := T, r (a,5,s) for convenience. Consider first equation (I5bj) 
with < k < K — 1. By the first variation formula and estimates of the size of the perturbation 
generating xj from S r (cr, 5, s), and calculating as in [TJ Prop. 27], we have 

ir 2 k+i{£ r (o-,5,s)) = / (H[El] - 2 - r 2 F(f r (o-,5,s)))xext,kJk 

= [ (7ro>^)- r2 / F(x,N Sk (x))))J k + 0(r*) 

Jdz«nsupp( Xext , k ) Xdx 1 Js k 
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= C 2 (sfc+i - e k ) - r 2 fi k {s,a) + 0(r 4 ) 

where v k is the unit normal vector field of dY$ n supp(xext,k) m S". 

Consider now equation (|5bp . In the neck we have if[E r (cr, 5, s)] = 0. Using similar estimates, 

ir 2k {£ r (a,8,s)) = - 2 - r 2 F(f r (a, S, s))) X neck,kh 

= "2 f Xneck,kh + 0(r 2+2 n 

= C^ef + 0(r 2 + 2 ») 



where 5 k is the displacement parameter of the k th neck. This is because l k is an odd function 
with respect to the neck having 5 k = 0, whereas the integral is being taken over the neck with 
5k i 1 0. Hence the integral f-^ nsupp ^ x k \ Xneck,klk picks up the displacement of the k th neck from 
its position at 8 k = 0. This same phenomenon arises in [TJ Prop. 27]. □ 

4.1 Proof of the Main Theorem 

It remains to find a value of the parameters (a, 5, s) so that £ r (a, 5, s) = 0. As shown in Lemma [H 
this is equivalent to find a solution of the equation B r (a, 5, s) = 0. In what follows, we will continue 
to assume that e = 0{r 2 ) and 5 = 0(r) and this will be justified shortly. As a consequence of 
Lemma [5l the equations that we must solve are as follows: 

C 2 £i = r 2 fjL 1 (a,s) + E[(a,6,s) 
Cx&x = Ei(a, 5, s) C 2 (e 2 - £i) = r 2 fi 2 (a, s) + E' 2 (a, 5, s) 

'■ and j 

Ck-iSk-i = Ek-i(<t,8,s) C 2 {e K -\ - ek-2) = r 2 fi K _i(a,s) + E' K _ 1 (a,5,s) 

-C 2 e K -i = r 2 fj, K (cr, s) + E' K (a, 5, s) 

where e k depends on a k in an invertible manner as indicated in Step 2 of the construction of the 
approximate solution, and E k , E' k are error quantities satisfying the bounds \E k \ = 0(r~ 1+2 ") and 
\E' k \ = C(r 4 ). We can abbreviate these equations by introducing the matrix M := (q j) where I 
is the {K — 1) x {K — 1) identity matrix and J is the K x (K — 1) matrix 



l 

-l l 



J 



-l l 
-l . 



The equations become 

M(8,ef = {E,r 2 n + E'f (6) 

where 5 := (5i, . . . ,5k-i), £ '■= (ei, • • • and so on for E, E' and fj,. 

We will solve these equations in two steps as follows. Note first that the matrix M is injective 
but not surjective, with vectors in the image of M satisfying the relation (0, e) • M(v, w) = for 
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all (v,w) G R 2K 1 , where e := (1, 1, . . . , 1). Let p : R 2K 1 — ► R 2K 2 be the orthogonal projection 
onto the image of M. The equation 

pM(e,5)=p(E,r 2 p + E') (7) 

can now be solved using the implicit function theorem. Thus the matrix pM : M. 2K ~ 2 — > M. 2K ~ 2 is 
invertible and the equation at r = is just pM{e, 5) = which has the solution (e, S) = 0. Hence 
the solution persists for small r. 

We now have a solution e := e r {s) and S r (s) of ([6]) for all sufficiently small r and depending 
implicitly on the one remaining free parameter s. Moreover, we see that e = 0(r 2 ) and 5 = 
0(r~ 1+2u ) = 0(r) since v E (1,2). It remains to solve ([7]) and we proceed as follows. Once (e,S) 
satisfy (J7|), then becomes equivalent to = (0, e) • M(e, 5) = r 2 e ■ p, + e • E' or simply 

^ MfcM*), s) + ^"(ffr(s), <5 r (s), s) = (8) 
fc=i 

where the error quantity satisfies the estimate \E"\ = 0(r 2 ). 

Equation ([8]) may or may not have a solution, depending on the nature of the function p^, 
which in turn depends on the specific nature of the prescribed mean curvature function F. However, 
if the following two conditions are met, then the implicit function theorem guarantees the existence 
of a solution. First, it must be the case that the equation at r = has a solution, in other words 
if the F-moments of the spheres Si, ... , Sk satisfy 

K 

k=l 

for some s, where p°(s) := (s + 2(k — 1),0, ... ,0). Second, if sq is the solution of this equation, 
then it must also be the case that the mapping 

K 

s^^MdBMi*))) 

k=l 

has non- vanishing derivative at s = so- If these conditions are satisfied, then the implicit function 
theorem implies that for r sufficiently small, there is a solution s(r) of ([5]). This completes the 
proof of the Main Theorem. □ 
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